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Abstract 



By studying the diffusion operator in Grover's quantum search algorithm, 
we find a mathematical analogy between quantum searching and classical, 
elastic two-body collisions. We exploit this analogy using a thought experi- 
ment involving multiple collisions between two-bodies to help illuminate why 
Grover's quantum search algorithm works. Related issues are discussed. 

PACS numbers: 03.67. 

Quantum computation is based on qubits. A qubit, literally a quantum 
bit of information, can be represented as a unit vector in a two-dimensional 
complex vector space. A typical vector can be written as \ip >= a\0 > 

>, with |ap + = 1. Measurement of a qubit yields jap and 
the probabilities of the system being in states |0 > and |1 >. Unlike a clas- 
sical bit, a qubit represents a superposition of the two states. This property 
leads to quantum parallel computation, which has been widely used in quan- 
tum algorithms [1]. Grover proposed a quantum search algorithm which can 
realize fast searching, such as finding an object in unsorted data consisting 
of items [2]. This algorithm can be described as follows. Initially, an n 
qubit system is set in an equal superposition of all basis states expressed as 



where the set {\i >} is orthonormal, and is usually far greater than 1. 
Each basis state corresponds to an item in the data. The problem is how 




(1) 



1 



to find the very basis state that corresponds to the object. This particular 
state is defined as the marked state, and the other states are defined as 
the collective state [3]. C and D are used to denote suitable inversion and 
diffusion operators, respectively. If C is applied to a superposition of states, it 
only inverts (i.e., changes the algebraic sign of) the amplitude in the marked 
state, and leaves the states in the collective state unaltered. D is defined 
a,s D = 2P — I, where P = (1/A^) Z^fjj-^o K J I- matrix notation, P 
is an N X N matrix whose entries are all and / is an x unit 

matrix. A compound operator is defined as U = DC. Each operation of U is 
called an iteration. After U is repeated 0{\/N) times, the amplitude in the 
marked state can approach 1. When a measurement is made, the probability 
of getting this state approaches 1 [2]. The algorithm can be generalized to 
the case where the marked state contains multiple basis states. [4] The search 
algorithm has become a hot topic in quantum information [5]- [7]. Grover 
outlined the steps that led to his algorithm [8], and proposed a classical analog 
of quantum search using a coupled pendulum model[9]. R. Josza provided 
another path toward understanding the nature of Grover 's algorithm[10]. In 
this paper, we discuss the similarity between quantum searching and classical 
two-body collisions and some related problems such as the limits on the 
algorithm. The similarity helps to understand why the algorithm works. 

We now examine some properties of D. If D is applied to a system in 
state Z^ilo^CjK >! it transforms the system into state Y^i^Jo^ di\i >. The 
relation between q and di is 

di = 2A- a, (2) 

where A is the average of all amphtudes of the system, namely, A = (l/N) Ylfjo^ Q 
Through easy computations, we find that 

N-l Af-1 

Y^di= Y^d (3) 

1=0 1=0 

and 

Af-1 Af-1 

y: ddi = E (4) 

i=0 1=0 

where c* denotes the complex conjugation of q. Eq.(3) shows that the trans- 
formation by D preserves the sum of all amplitudes. Eq.(4) results from 
normalization of the wave function. 
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We suppose that the collective state contains A^i basis states with equal 
amplitude Qq and the marked state contains N2 basis states with equal ampli- 
tude 60 before U is applied. From now on, we assume that all amplitudes are 
real numbers. C transforms Bq to —Bq, and leaves oq unaltered; D transforms 
—60 to bi, and ao to ai. According to Eqs. (3) and (4), two equations 

A^iai + N2bi = Niao - N2bo, (5) 

Nial + N2b\ = Nial + iVs^o (6) 

are obtained. They remind us that the operation of D is analogous to the 
elastic collision of two bodies. We assume they are two rigid balls with masses 
mi and Ball 1 and ball 2 have velocities u and v, respectively. Before 
the first collision, these velocities are indicated by Uq and —Vq; after the 
collision, they are denoted by ui and vi. All velocities are confined to one 
straight line and the minus sign means that the two balls move in opposite 
direction before the collision. We obtain two equations 

niiUi + 777,2^1 = miUo - m2fo, (7) 



1111 

2,-'- 2-"- 2,-'- 2 /o\ 

based on the conservation of momentum and mechanical energy [11]. We 
find the forms of Eqs. (7)- (8) similar to the forms of Eqs. (5)- (6). Their 
solutions also have similar forms. The solutions for Eqs. (7) and (8) are 

(mi - m2)Mo - 2m2Vo 
ui = ■ , 9) 

mi + m2 

2miMo + (mi - m2)vo 
vi = ■ . (10) 

mi + m2 

Let Vc = {rriiUQ — m2f 0) / {frii + m2) be the velocity of the mass center of the 
system. Vc is constant during the collision of the two balls, ui and vi can 
thus be represented as 

Ml = 2Vc - Uq (11) 



vi = 2vc + vo (12) 
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Now, we assume that ball 1 and ball 2 consist of A^i and N2 particles 
respectively. For convenience, we denote these particles aspiii = 0,1, N — 
1). The mass of each pi is denoted by = mo, where mo denotes unit 
mass. Therefore mi = Nim^, m2 = N2mo, and iVi + A^2 = N. Because of the 
similarity between Eqs. (7)-(8) and Eqs. (5)-(6), a one-to-one correspondence 
exists between the amplitude of |z > and the velocity of Pi. All velocities are 
confined to one straight line because all amplitudes are assumed to be real 
numbers. The probability of getting |i > is analogous to the kinetic energy 
of Pi. Ni particles in ball 1 and N2 particles in ball 2 are analogous to Ni 
basis states in the collective state and N2 basis states in the marked state, 
respectively. Because Vc is also represented as Vc = (l/N) Y.f=o^ where Vi 
is the velocity of Pi, it can be concluded that Vc is analogous to A in Eq.(2). 
We summarize the above-mentioned correspondences in Table 1. 

Table 1: Correspondences between Grover's search algorithm and classical 

two-body collisions 



Grover's search algorithm 


classical two-body collisions 


Eqs. (5)-(6) 


Eqs. (7)-(8) 


state \i > 


particle pi 


amplitude of i > 


velocity of Pi 


probability of getting \i > 


kinetic energy of Pi 


marked state 


ball 2 


collective state 


ball 1 


A ( average of amplitudes) 


Vc (center of mass velocity ) 



Initially, two balls move with the same velocity, i.e., all particles move 
with the same velocity. This is analogous to the initial equal superposition 
of states in quantum searching. The operation of C is analogous to the 
operation that inverts the velocities of particles in ball 2, and leaves 
particles in ball 1 unaltered. From Eqs. (11) and (12), in the center of mass 
frame, we find that before the collision, two balls have velocities uq — Vc and 
—vq — Vc] after the collision, their velocities are transformed to Vc — uq and 
Vc + vq- This fact means that the velocities are inverted by the collision. In 
the reference frame of the laboratory, however, the after-collision velocities 
are 2vc — uq and 2vc + vq, respectively. Application of the diffusion operator 
D, which is also called inversion about average, transforms the amplitudes 
of the collective state and the marked state from ao and —bo to 2A — and 
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2A + Bq. Therefore, the similarity between operation of D and the colhsion 
is obvious. 

According to the classical analog of the quantum searching algorithm, we 
conceive a thought experiment to simulate the algorithm. Ball 1 and ball 2 
move rightwards along a horizontal and straight path without friction with 
the same initial velocity uq = vq = v as shown in Fig. la. In order to invert 



Uo=V Vo=V V -V Ui Vi 




Figure 1: A scheme used to simulate the quantum search 
algorithm using classical collisions The two balls of masses mi 
and 777,2 are moving with the same velocity along a horizontal and 
straight path without friction (Fig. la). After ball 2 collides with 
the obstacle at the right end of the path, its velocity is inverted 
(Fig. lb). After the collision of the two balls, the velocities become 
Ui and Vi (Fig.lc). 

the velocity, we set an obstacle at the right end of the path. The velocity 
of ball 2 is inverted but its modulus remains unaltered after it collides with 
the obstacle. The velocities of the two balls after the collision are v and —v 
as shown in Fig. lb. This collision is analogous to the operation of C used 
in the quantum searching algorithm. After ball 2 collides with the obstacle, 
it collides with ball 1. The second collision is analogous to operation by D. 
The velocities after this collision are Ui and Vi as shown in Fig.lc. The two 
collisions that are analogous to C and D constitute the complete iteration 
analogous to U used in the quantum searching. There are 3 possible cases 
after the collision of the two balls: 1) The two balls both move rightwards, 2) 
ball 1 moves leftwards and ball 2 moves rightwards, and 3) the two balls both 
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move leftwards. In some cases, the next collision of two balls cannot happen 
after the velocity of ball 2 is inverted again. Because we only are concerned 
with collisions and have no interest in the positions of the two balls, we can 
exchange the two balls' positions or set an obstacle at the left end of the 
path in order to make the iteration continue. Assuming the two balls have 
velocities Un and w„ after n iterations, we obtain the recursion equations in 
matrix notation 

= [ ^ r I (13) 

V N1+N2 N1+N2 / \ ^ / 

by replacing the subscripts 1 and in Eqs.(9) and (10) by n + 1 and n, 
respectively, and by using the conditions that mi = Nirrio, and m2 = A^2^o- 
For the quantum system, after U is repeated n times, the amplitudes in the 
collective state and marked state are transformed to a„ and 6„. One can 
directly obtain the recursion equations 





N1-N2 -2N2 
N1+N2 

2Ni 



N1+N2 N1+N2 




(14) 



from Eq.(13), based on the above-mentioned similarities. The initial condi- 
tions are Uq = Vq = v for Eq.(13), and oq = &o = ^/VN for Eq.(14). The 
2x2 matrix in Eq.(14) is the product of Grover's diffusion operator and 
suitable inversion operator, and is denoted by 



jVi-jV2 -'2N2 
2Ni 



^1+-^^ ^+^, . (15) 



N1+N2 N1+N2 



We will proceed to solve Eqs.(13) and (14) using linear algebra [12]. One 
finds that 

(16) 



Un 


j = T" 1 








[ Vq 


an ^ 


j =T" 1 


( ao 


K , 


[ bo 



(17) 

where n = 1, 2, ■ ■ -. In order to obtain explicit expressions for m„, Vn, an and 
hni we transform T to the diagonal matrix 

Tj, = S~'TS ^ f V M • (18) 
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The eigenvalues X± of the matrix T are the solutions of det(T — XI) = 0. 
They can be expressed as 

\ _ Nl-N2±2iy'NlN2 (-IQ\ 

N1+N2 ■ 

5* and 5*"^ are expressed as 

^= ( . . k V (20) 

^'^ = ^1 : ~'7k I- (21) 



One can easily find that 

T" = STJ^S-\ (22) 
If = 1, and N ^ 1, Td can be expressed as 

Tn = ( '7 ) . (23) 

using tan 2^ ^ 26, where 9 = Using Eqs.(20)-(23), one can derive that 

rn^f cos(2n^) -^sin(2n^)\ 

V yiV^sin(2n^) cos(2n^) / ' ^ ^ 

Using Eqs.(16), (17) and (24), we obtain the explicit expressions 

%cos{{2n + m \ ^25) 



Un \ V- 



Vn J \ vy/Nsin{{2n + l)9) 

( <^n\^( :;hi^os{{2n + m\ 

[bnj [ sinii2n + l)e) J' ^ ^ 

noting that sin 6* ~ and cos 6' ~ \f^^- During the course of colli- 
sions, energy is transferred between the two balls through the change of their 
velocities. Let riQ «i (tt\/N/A — 1/2) be the integer obtained by rounding 
(tt^/N /A — 1/2). When n = tiq, ball 2 can acquire almost all the mechanical 
energy of the system. Correspondingly, for the quantum searching algorithm, 
the amplitude in the marked state can approach 1 if the number of repetitions 
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of U is Uq. The quantum system lies in the marked state if a measurement 
is made. Our results show agreements with Refs.[4][13][14]. We now discuss 
some related problems based on these discussed similarities. 

1. Simulating Grover quantum searching algorithm on a classical 
system Because pi is analogous to \i >, a system consisting of 2" particles 
Pi{i = 0,1,2,..., 2" — 1) can be used to simulate a quantum system that 
stores n qubits of information. For example, if mi = 3mo and m2 = rriQ, the 
classical system can simulate a two qubit system with A'^i = 3, = 1, and 
= 4. After exactly one iteration, we find that ui = 0, vi = 2v. This result 
means that ball 2 obtains all the energy of the system, which is analogous 
to the probability equal to 1 of getting the marked state. Our result is in 
agreement with the result demonstrated on an NMR quantum computer [5]. 
If mi = 7mo, and m2 = mo, the following results are obtained: 



Ball 2 nearly gets the entire energy of the system after 2 iterations. Because 
mass is continuous, the case where m2 = mo, and mi = NirriQ is equivalent 
to the case where m2 = kmo, and mi = kNimo {k is a positive integer). The 
case where mi = 3mo, and m2 = mo is analogous to the case where A^2 = N/A 
as discussed in Ref . [4] , in which after one iteration, the probability of getting 
the marked state is 1. If we choose the condition that mi/m2 = p/q, where 
g 7^ 1, and p and q are positive integers with no common factor, we can 
simulate the quantum search algorithm in the case where the marked state 
contains multiple basis states, which cannot be realized by NMR [15]. 

2. Limits on quantum searching Because the two balls exchange 
their energy after they collide with each other, neither of them can get all 
the energy of the system if mi = m2. This case is analogous to the case 
where the search algorithm is invalid, more specifically, where Ni = N2 [16]. 
We also find that 



using the initial condition uq = Vo = v. Eq.(29) means that if N2 > N/4, the 
velocity of ball 1 is inverted after the first iteration. It is possible that the 




(27) 




(28) 




(29) 
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two balls move in the same direction after the velocity of ball 2 is inverted 
in the second iteration. The energy of ball 2 is reduced after it catches up 
with ball 1 and collides with it. In this case, the searching algorithm is not 
efficient [16]. 

3. Qualitative discussion of the needed number of iterations Con- 
sidering the case where m2 = rriQ, mi = {N — l)mQ. The initial condition 
is Uq = Vq = V, and ^ 1. We obtain ui ^ v, and Vi = 3v after the 
first iteration. This means that the velocity of ball 2 increases by 2v and 
the velocity of ball 1 hardly change in any iteration. The velocity of ball 2 
increases by 2nv after n iterations. When n ~ /2, ball 2 gets almost all 
energy of the system. 
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